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The phenomenon described by our title should surprise no one. What may be 
surprising though is how easy it is to produce a quantum system with this feature; 
moreover, that system is one that is often used for the purpose of showing how sys- 
tems equilibrate. The violation can be variously manifested. In our detailed example, 
bringing a detuned 2-level system into contact with a monochromatic reservoir does 
not cause it to relax to the reservoir temperature; rather, the system acquires the 
reservoir's level-occupation-ratio. 



I. INTRODUCTION 

The zeroth law of thermodynamics is a kind of transitivity. Suppose that if system A is in 
contact with system B they do not exchange energy. Similarly if A is in contact with another 
system C, they also do not exchange energy. Then if B and C are brought into contact they 
too will not exchange energy. As suggested by its position among the laws, this observation 
about Nature is what permits the other laws of thermodynamics to be discussed. 

In this article we present a system that violates the zeroth law. It is a model, and the 
actual physical systems for which this is (sometimes) an excellent approximation would 
presumably satisfy the zeroth law on time scales for which the modeling broke down. Nev- 
ertheless, there may be experimental time scales for which an apparent violation could be 
exhibited. In any case, like Arnold's cat map jlj, this not-quite-physical system serves to 
elucidate the foundations of thermodynamics. 

The key to the breakdown is the fact that our model system is not ergodic. In classical 
mechanics ergodicity is usually phrased as the equality of time and phase-space averages. 
This is equivalent to the non-existence of any constant of the motion except energy (we 
assume the system is confined so as to fix momentum and angular momentum). In the 
quantum context, we take ergodicity to mean the absence of non-energy constants, and 
what lies behind the phenomenon we exhibit is just such a quantity. 

A fundamental property of temperature (see for example [2j) is that objects placed in 
contact are not in equilibrium unless they have the same temperature. In demonstrating 
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the zeroth law breakdown, we will give an example in which two systems at the same 
temperature, when placed in contact, evolve to different temperatures. Call one system 
"A" (which will be thought of as a reservoir) and the other "S," a two- level system. After 
contact (as described below) they will be at different temperatures. If A is then placed in 
contact with a second system (say S') that has energy spacing different from that of S, S' 
will evolve to a third temperature. However, if S and S' are brought in contact, they will 
exchange energy, violating the zeroth law. 

For the purpose showing the approach to the non-equilibrium stationary state we will 
use an intermittent contact. However, having attained that state, reservoir contact becomes 
arbitrary; it can continue to follow our scheme or it can remain indefinitely coupled, with the 
density matrix of the system remaining fixed and out of equilibrium. To illustrate subtleties 
of our contact process we will treat two interaction schemes. For both of them the systems 
evolve to temperatures different from the reservoir, and for the one that does not violate 
the zeroth law the specifics of the contact indeed play a role. 



II. MODELS AND NUMERICAL RESULTS 

As indicated we work with two interaction models, one for didactic purposes only. The 
two interactions are commonly called "rotating wave," which will exhibit the zeroth law 
violation, and "spin-boson," which as stated, is included for purposes of contrast. The 
Hamiltonians on which we focus are 

H = H A +H S + 1 V 1 H A =Y,^a\a kl H s = EsU, (1) 

k 

where 

V = [gkO-kS* + g*k a { s ) ) rotating wave (RW), (2) 

k 

or 

V = (V + s) gk + QfeJ , spin-boson (SB) . (3) 

k 

Our notation is standard. System A is a collection of harmonic oscillators (thought of as 
photons or phonons) and a) and a are its creation and annihilation operators. For system S 
we have similar operators, and s, although since we take S to have only two levels, squares 
of these operators are zero. The quantity 7 multiplying the potential is a convenient overall 
scaling. Both interactions are physically significant. For RW, using the Wigner-Weisskopf 
method it is shown in that if the oscillator states have a thermal distribution with inverse 
temperature [If, then the ratio of probabilities for the excited and ground states of S is 
e~^ E — the usual equilibrium. RW dynamics has, besides H, the constant of motion 

N = a[ak + s^s (excitation number) , (4) 

k 

and is not ergodic. From the demonstration 0] it would seem that there are no thermody- 
namic consequences. 

This satisfying situation is the result of approximations and assumptions on A. Consider 
a single frequency reservoir that is detuned from S, i.e., E 7^ 00. We want A to act as 
a reservoir. It begins at temperature 1/(3 and couples to S for a short time, after which 
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they are decoupled; A is refreshed and they are re-coupled. Iterate. This models the 
Markovian assumption in studies of dissipative dynamics in that A loses its correlations 
with S. Similarly, in diffusion a particle connects and disconnects in successive encounters 
with other particles (viewed as a reservoir), with correlations carried away by translational 
degrees of freedom. Our coupling-decoupling has the following mathematical realization: let 
the initial density matrix be p(0) = ps(0) <g) pa(0) with pa(0) = (l — e _/3w ) exp (—/3Ha)- 
Evolving under H, p(At) = f/p(0)t/^, with U = exp(-iHAt). Next separate the systems 
and form the reduced density matrix for S, ps{At) = Tr^p(At). The nature of A as a 
reservoir is manifested by assuming that it recovers from its interaction with S, and is again 
in equilibrium at temperature-l//3. Repeat this many times. Thus the time-2At state of S is 
given by ps(2At) = Tr aU {ps(At) <g> Pa(0)} U'. A physical realization could be the passage 
S through a sequence of sharply tuned cavities — but not tuned to the frequency of S. The 
coupling is weak, but one expects S to eventually acquire the temperature of A; it does not. 
We find, both numerically and analytically, that the relative probability of being excited is 
exp(— {3cu), not exp(—f3E). 

Remark: We speak of S's disequilibrium in terms of temperature for convenience only. For a 3- 
(or more) level system there is also a failure to reach the appropriate Boltzmann distribution. 
Moreover, if one uses the rotating wave coupling, then the zeroth law breakdown described 
in this article does obtain. In particular, S adopts the level-occupation probabilities of 
the reservoir and it takes no energy to maintain this disequilibrium. However, only if S's 
levels are equally spaced can one still speak of temperature. (A similar result in different 
circumstances was found in (5)].) 

The disequilibrium, however surprising, is not by itself a consequence of non-ergodicity. 
The SB model, without the constancy of N, behaves similarly, although the final temperature 
depends on other variables. The intermittent interaction means that the overall effective 
Hamiltonian is time-dependent, hence need not conserve energy. Thus when contact is 
broken there can be interaction energy, ^yTr(p(At)V), whose elimination requires external 
work and thus maintains the disequilibrium [6|. Where the RW model is different is that— 
as we will show — in the non-equilibrium stationary state there is no A-S decoupling energy. 
Separating the systems costs nothing and nevertheless disequilibrium is maintained. Since 
energy is conserved during the contact interval, the coupling energy equals (up to a sign) 
the change in (Ha + H s ). Thus if S has reached its steady state, the cost of maintaining 
disequilibrium is the change in (Ha) ■ For RW this cost is zero. After quantitative analyses, 
we will give a qualitative argument for RWs peculiarity, with the role of non-ergodicity 
evident. 

In Fig. [TJ is shown the time-dependence of /3-effective, computed from the occupancy 
ratio, for a 2-state RW system for which E = oj/2. The initial /3 for S is 1, so that as 
they repeatedly touch and pull away, S cools. Fig. [2] shows the interaction energy. As S 
approaches its (non-equilibrium) steady state, the RW interaction energy goes to zero [7|. 

In the sequel we first present a density matrix formalism in which the standard equilibra- 
tion result is shown to be approximate. Then we specialize to single-frequency reservoirs, 
where the evidence is unambiguous: For RW, S adopts the occupation ratio of the bath, not 
its temperature. If A has frequency cu, that ratio is exp(— f3u), giving S an effective tem- 
perature T e fj = (£'/a;)Tbath- This is true whether the A-S contact duration is long or short 
and whether or not S begins in a diagonal density matrix. For SB, S tends to an effective 
temperature that is (in general) not that of the bath. For a certain limit of 7 and At the 
temperature of S tends to infinity. Moreover, for RW maintaining disequilibrium costs no 
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FIG. 1: Effective-/? on successive contacts with the (3 = 1 reservoir, for the RW interaction. Initial 
conditions as well as duration and strength of the interaction vary. For Run-2 ps(0) is not diagonal. 
In all cases the same non-equilibrium limit is attained. 
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FIG. 2: Interaction energy at the end of the A-S contact period for the RW interaction. Coupling 
parameters and initial conditions vary. For Run-2 ps(0) is not diagonal. 

work, i.e., the coupling and decoupling requires no energy once S has become stationary. 
However, for SB, in its nonequilibrium steady state, during each decoupling the bath gains 
energy, an increase that comes at the expense of whatever external process implements the 
decoupling (cf. [jjj]). 

III. SMALL TIME CALCULATION 

For a general reservoir, single or multiple frequency, we begin with p(0) = ps(0) <S> Pa, 
where ps(0) is arbitrary and pA = e~^ UA /Za-, Za = Tre _/3H ^. For small At 

At 2 

p(At) = p(0) - iAt[H, P (0)] - —[H, [H,p(0)]\ + 0(At 3 ) . (5) 

In analyzing Eq. ()5]) we use the following: 1) [Ha,Ps(^) ® Pa] = 0; 2) the diagonal matrix 
elements of V on the eigenvectors of Ha are zero (true for both RW and SB). After several 
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steps Eq. fl5]) reduces to 



Ps(At) = p s (0)-iAt[H s ,p s (0)} 
Af2 [H s , [Hs,p s (0)}} 



2 

7 2 At 2 



Tr A [V,lV,p s (0)®p A }] + O(At 3 ). 



(6) 



Although our results are more general, we specialize Eq. OH]) to a "diffusion limit" [9| , namely, 
7—7-00 and At — > such that 

7 = 7V / At — > const . (7) 
In this limit we compute [p^(At) — Pa(0)] /At, leading to 



dps 
dt 



-i[H s , ps{0)) - ^ Tr A [y, [V, p s ® Pa] 



We first apply these results to RW. Calculating the appropriate matrix elements and 
commutators, we obtain two independent equations. For the diagonal 



dp 



s,oo 



dt 



k 



\9k\ 



[e^Psfio - Ps,u] 



(8) 



The off-diagonal element satisfies 

dp s< w _ 
dt ~ 



^2 



iE-i-y 
2 ^ 



9k\ 



1 + e~P" k 



PS, 10 



(9) 



(and p 5i0 i = ps.io*)- As a result \p s . 



01 



r = yJZ Is-fcl 2 coth 



e (iE r)t q ? w h ere 



Using ps n = I — p5 j0 o m Eq. flE]) yields for the limit 

" s '°° (oo) - £ fc M 2 coth(fc) ■ 

If g% is sharply peaked at ko such that E = Uk , a short calculation shows that 



(10) 



Ps,n(oo) 
Ps,oo{oo) 



-0E 



(11) 



This recovers thermal equilibrium, and is the result of Einstein [10] • However, this conclusion 
depends on the {gk}- In general, one does not recover thermal equilibrium. Fig.[3]shows how 
effective temperature varies with the width of a Gaussian approximation to a 5-function. 
The temperature is off, albeit not grossly. 
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FIG. 3: Variation of temperature as the spread of a Gaussian pseudo-5-function is varied. Tem- 
perature, T, for S is computed from (1 — ps,oo)/ Ps,oo = exp(— E/T) with ps.oo given in Eq. (fl"Qj) . 
(For a Lorentzian pseudo-rj- function the deviation is larger.) 



For the monochromatic, detuned thermal bath, "g^" is a ^-function at some k, but now 
E 7^ ujk- Eq. ({TO]) yields the probability ratio e~ l3uJk , the dropoff exponential of the bath. As 
a result S has effective temperature 

E 



-off 



r, 



bath " 



or (3 eS = (3 



E 



(12) 



We turn to the SB interaction. Calculating the matrix elements in the limit ([7j), one 
obtains 



dps, 



oo 



dt 
dps,w 



I \9k\ 2 coth 

\ k 



(3uk \ 



2 J [ps,oo - Ps,n] 



dt 



iEpsfii-1 > Jflfc| coth — - [ps-oi - p s ,io 



It follows that p5 o — > \ and p^oi — )• for t — > oo. The system tends to infinite temperature, 
which is consistent with [8(, although the recoupling scheme differs. 



IV. ROTATING WAVE, ARBITRARY TIMES 

Returning to RW, we show that Eq. (jT2l) holds for any interaction duration. The single 
frequency reservoir only has transitions within "H n = {|n, 0), |n — 1, 1)} (1 st label A, 2 nd S), 
for n > 1, a manifestation of non-ergodicity. Taking g real (for simplicity) and 7 = 1, the 
Hamiltonian is 

nu gy/n 
g^/n (n — l)u + E 

= (nco + ^)l + a z (^) + (13) 

with the ex's the Pauli spin matrices. Defining Vl n = y ((u — E)/2) 2 + ng 2 , rh n = 
(g^/n, 0, (u; — E)/2) /O n , and f2 = nu + (a; + E)/2, the propagator in is 

[1 cos Q n At — ier ■ rn n sin Q n At] . (14) 



e -iH„At _ e ~in At 
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Assume that ps(0) is diagonal in the energy basis, so that transitions away from each \ua, n$) 
can be analyzed separately (the general case is treated in Appendix |A|) . From Eq. (jl~4"]) , if one 
begins in S's ground state, the probability of excitation is r n = {ng 2 /O 2 ) sin 2 Q n At (< 1). 
For de-excitation the probability is the same, as is evident from Eq. (1141) . Calling the ground 
state probability, summed over all ua, q, the change in q after a single interaction interval is 

= -q^2p n r n + (1 - q) J^Pn?Vt-i • (15) 

In the second sum r takes the index (n + 1) because these contributions to q come from 
% n +i- Their weight in pA is nevertheless p n . Both sums can be taken from to oo since 
r = 0. The second sum can be written ^2p n r n+ i = J2Pn+i r n+i- Define T = J2Pn r n- 
Then Eq. ffT5l) becomes 

Aq = -qT + (1 - q)e~ puJ T . (16) 

Setting Aq = to find the stationary state gives (1 — q)/q — e _/3a \ (It is also easy to see 
that convergence to this value is exponentially rapid.) 

The same argument can be applied to changes in A's probability distribution during the 
course of S's coupling. The possible transitions and the probability transfer for each are 
as follows: Losses in p n : [\n, 1) — >• \n + 1,0), r n+1 ], [\n, 0) — > \n — 1,1), r n ]. Gains in p n : 
\n + 1,0) -)■ \n, 1), r n+ i], [\n - 1, 1) -)■ |n, 0), r n ]. Thus Ap n = - [p n (l - g)r n+i +p„gr n ] + 
[p n+ igr n+1 +p n _i(l — q)r n ]. The limiting value of q is given by (1 — q)/q — e~^ w . Using 
the fact that initially p n +i/p n = e _/3a; , we find that Ap n = 0. When S has reached its non- 
equilibrium steady state there is no change in the reservoir's probability distribution, hence 
no change in its energy. Moreover, by conservation of energy during the period of contact, 
there is also no residual coupling energy. (This is also true for a non-diagonal ps(0). See 
Appendix |A~1) 

V. SPIN-BOSON: DISEQUILIBRIUM TAKES WORK. 

For SB, numerically, the temperature does converge, but to a value that depends on 
details (see Fig. H]). In the diffusion limit, the propagator becomes U = exp(— i^yVy/At). 
The SB interaction induces the following transitions: for n > 0, (n, e) — > (n + 1, 1 — e), 
e = or 1, with matrix element 'yy/Aty/n + 1, and for n > 1, (n, e) — > (n — 1,1 — e) 
with matrix element ^\f~Aiy/n. As in Eq. f|T5|) . we monitor changes in q. Now 
both sums have the same range and using the squares of the matrix elements, Aq = 
■y 2 At{—qJ2Pn [n + (n + 1)] + (1 — q) J^Pn [n + (n + 1)]}. Requiring Aq = implies q = 
1 — q, or equal probability of excited and ground states, namely infinite temperature. How- 
ever, unlike the RW case, the reservoir does not go to its original state at the end of each 
connection interval. We have calculated, both analytically and numerically, that the energy 
cost in this limit is 7 2 At, so that per unit physical time there is a cost, 7 2 , to maintain the 
disequilibrium. 

VI. DISCUSSION 

For RW, we provided a simple argument for the equality of A and S occupation ratios. 
But even in our "simple" argument it may not have been clear where non-ergodicity entered. 
In the following we do it in words. When S is in its steady state each contact with A must 
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time 

FIG. 4: Effective-/? on successive contacts with the (3 = 1 reservoir. Duration and coupling strength 
of the SB interaction vary. For the curve using "x" symbols ps{0) is not diagonal. The two lowest 
curves are for decreasing times and fixed 7, showing the approach to /3 e ffective = 0. Numerically, 
for small At, (3^ cAt 2 , with c 0.15. 

leave it with as much energy as it had before. This does not imply (cf. SB) that the final 
energy of A is unchanged. However, with the additional conservation law of total excitation 
number (Eq. (jl])) for RW, the reservoir must also be left with its original excitation. Energy 
constancy for the system implies constancy for the reservoir. (This argument also applies if 
S has more than two levels.) 

With a multi-frequency reservoir this argument fails, since energy transfer can be ac- 
complished by superpositions of excitations of many frequencies. Moreover, for any finite 
number of reservoir oscillators greater than one, the decoupling process does take energy 
(see Appendix IB"]) . 

For both models the systems are in contact with a monochromatic reservoir and have tem- 
peratures different from that of the reservoir. By itself this violates nothing: the connection- 
disconnection process allows a transfer of energy that can maintain disequilibrium. Where 
the surprise lurks is that for the rotating wave model, this transfer is zero. The systems 
are in intermittent contact, but S is in a state of disequilibrium that takes no energy to 
maintain. Note too that once stationarity has been attained the contact can be intermittent 
or not, since our effect is independent of At. Another way of phrasing the violation is that 
after each contact period it could happen (and does happen for SB) that the final reservoir 
state has changed its energy in such a way as to maintain the disequilibrium. However, for 
RW this "change" is zero. 

In principle, the effect we have reported here must disappear on a sufficiently long time 
scale, since in the real world the rotating wave interaction is only an approximation. How- 
ever, if the detuning is slight it will be a good approximation, so that a variation of the 
physical setup mentioned earlier (an effectively two-state system passing through a succes- 
sion of detuned cavities) could yield the disequilibrium discussed here. 
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Appendix A: Exact computation of density matrix evolution, rotating wave (RW) 

model 

Tin is the subspace spanned by {\n, 0), \n — 1, 1)}, n > 1, 1 st label A, 2 nd label S. Because 
of the non-ergodicity of RW dynamics, transitions can only occur within a given % n and 
one can consider the following restricted Hamiltonian (taking g real and 7 = 1): 



H„ = 



g\fn in — l)u + E 



( u + E\_ t (u-E\ r - 
= I noj H — I 1 + a x I — - — J + gVn<r x , (Al) 

Define Q n = \J {{00 - E)/2) 2 + ng 2 , rh n = (g^/E, 0, (u - E)/2) /ft n , and Q = nw + 
The propagator within T-L n is 

U n {t) = exp (—iH n t) = exp(-iQot) [1 cos Q n t — ier ■ m n sin Q n t] 
cosf2„t — sinf2 n t —i^^smil n t 



— i%^sinfi n i cosf2 n t + i^^^- sin il n t 



The last line implicitly defines c n , etc. 

The general form of the initial density matrix is 



P(°) = $Z Pn l^fi + P) n n,i] + Pod , (A3) 



n>0 



where p n = (1-e #")e /3nw , 7r raAi „ s = |n4)(n A ||ns)(ns|, n A = 0, 1, 2, . . . , n s = 0, 1. p is real 
and < /i < 1. The term poo is the off-diagonal (for S) contribution. It is given by 



PoD(0) = {C\0s)(ls\+C\ls)(0s\}Y l P"\ n )( n \> with|C| 2 <i-(p-i) . 

n>0 "V / 



(A4) 



The condition on ( is guarantees that p(0) is a bona fide density matrix. 

The time-dependence of p is obtained from p(t) = Up^W. This can be evaluated for 
each term in the sum and gives the following for p(t): 

pit) = p pn 00 + ^2p n ^U n TT nfi Ul + - n)U n+1 n n:1 U n+1 + U Pod{0)U^ . (A5) 

n>l n>0 
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where U n is the matrix given in Eq. (1A2j) (or its equivalent in bra-ket notation). In Eq. (1A5j) 
U n+ i appears in the second sum, since this acts in "H n+ i. Note too that % Q (with basis |0, 0)) 
is constant in time. The off diagonal term is 

Pooit) = Up OD (0)rf 



n 00 + 2j (^,0 + 7Tm-l,l) U ] m 
m>l 



j ' v° v ' v° °y ' ^ 1 

We evaluate the time evolution. First the projection operators: 



X>HHC|o s )(i 

+ adjoint , (A6) 
This becomes 

PoD(t) = ( 1 ^2pnU n \n,0)(n,l\Ui +1 + adjoint. (A7) 

n>0 

In each subspace H n , S's operators can be written as 2-by-2 matrices 

|0><0| = (I ° n ), \1)(1\= (° n ?V |0><1|= (° n l), |1)(0|= ( ? °Y (A8) 



= ^n^o + (1 - ZnHn-i.i + (^„-R„ + adjoint) , (A9) 
with R n = \n, 0)(n — 1, 1|, the raising (with respect to A) operator. Similarly 

n jj\ = ( 7n+l s n+l -Z7 n+ iS n+ i(c n+ i - z5 n+ iS n+ i) 

n+1 " +1 " \hn+lSn+l(c n+1 +i5 n+ iS n+1 ) C 2 n+1 + 5l +1 S 2 „ +1 

1 ~~ ^ra+l —Zn+1 



Z n+l x n+l 

= (1 - a;„ + i)7r n+ i i0 + x„ + i7r nil - {z n+1 R n+1 + adjoint) . (A10) 
The off-diagonal term mixes % n 's. We look at 

PoD{t) = Cy^Pn^nK Q) (fl"n+iK l)) 1 " + adjoint . (All) 

n>0 

We write this in full: 

Pooit) = C^Pn[C/n(l,l)K0) + C/ n (2,l)|n-l,l)][?7 n+1 (l,2)|n + 1,0) + C/„ +1 (2,2)|n,l)] f 

n>0 

+ adjoint 

= C^Pn{^n(l,l)|n,0)(n+l,0|t/ n+1 (l,2)* + C/ n (l,l)K0)(n,l|C/„ +1 (2,2)* 

n>0 

+U n (2, l)\n - 1, l)(n+ 1, 0|[/ n+ i(l, 2)* + C/ n (2, l)|n - 1, 1) (n, l\U n+1 (2, 2)*} 
+adjoint . (A12) 
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When performing Tr^, terms with different tia values go away. Substituting the relevant 
quantities from the expression for U n one gets 

)\n, 0)( iS n+1 s n+1 )] + adjoint 

n>0 

+ [terms that disappear under Tr^] . (A13) 
Substituting all this into Eq. (IA50 yields 

p(t) = p prr 00 + ^p n p {x n 7i nfi + (1 - x n )7r n _i i i + {z n R n + adjoint)} 

n>l 

+ ^2Pn(l ~ p) {(1 - x n+1 )n n+lfi + x n+ i-K n .i - (z n+ iR n+1 + adjoint)} 

n>0 

+p OD (t) . (A14) 
with Pod(^) given in Eq. (IA13I) . Changing dummy indices on some of the sums leads to 

+ ^Pn+iP {z n+1 R n+ i + adjoint) 

n>l n>0 n>0 

+ ^2Pn-l(l - /i)(l - X n )n nfi + 2jp n (l - p)x n+1 7T ntl 

n>0 

^Pn(l - P) {z n+1 R n+1 + adjoint) + PodOO • (A15) 



n>l 



n>0 



Our primary interest is Ps(t) so we consider what happens when Tr^ is applied: 



Ps(t) = |0)(0| 



PoP 



+ ^PnPXn + J~]p n -lQ- -//)(!- X n ) 



n>l 



n>l 



^p n+ ip(l - x n+ i) + 2jp n (l - p)x n+1 

_n>0 n>0 

+ 1 ^5Z Pn t( c « -^n s n.)l°)( 1 l( c n+i + adjoint I . (A16) 

I n>0 J 

Now although mathematically c n and c n+i can simultaneously equal one, physically it is 
not reasonable, so that generically |(c n — iS n s n )(c n+ i — iS n+ is n+1 )\ < 1. This implies that 
the off-diagonal terms drop away. We also use the fact that p n+ i = e~^ u p n . Without the 
off-diagonal terms, Eq. (1A16j) becomes 



Ps(t) = |0)(0| 



POP + ^Pn [pX n + e^(l - p)(l 



n>l 



n>l 



(A17) 



12 



We check the Aq calculation of the main text. Eq. (1A17|) gives the new value of "q" as 



q = Poq + ^Pn [qx n + e^(l - g)(l - x n j\ (A18) 



n>l 



Subtracting q = q J2 n >o Vn from both sides gives 

Aq = {po - Po )q + J2 Pn W x « ~ X ) + ~ ^X 1 ~ 2 

n>l 

= -5^(1 - *n) [q - e*"(l - q)} , (A19) 



n>l 



implying that e & u is the excited to ground state occupation ratio in the stationary state. 

Appendix B: Energy flow in the diffusion limit 
1. Energy flow in S 

At time t, we call 

E s (t) = Tr (p s (t)H s ) = Ep s ,u(t) • (Bl) 
Then in the diffusion limit the flow of energy of S is 

dE S _ OT T7„-2Tt 



dt 



2TEe- 2Lt [ps,iM-Ps,n(oo)} (B2) 



which goes to zero. 



2. Flow of energy in A 



At time t one has 



dE A l_. 



tjf Tr (H A [V, [V, p s (t) ® p A \\) (B3) 

where p A is the thermal state of the bath, or after explicit calculation when t — > +oo 

dEA -|7 2 V \gk\V Hk) - u(£)) {e~^ - e~^) 

dt Efck|2coth ^2. (l- e -M*))(l- e -^M) 1 } 

In particular, > for t — > +oo unless the bath has only one frequency (tuned or 
detuned). In that case ^ ->■ 0. 

Therefore for a multi-frequency bath A, the bath will absorb energy at a constant rate, 
the energy coming from whatever external force is inducing the intermittent contact. 
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3. Flow of interaction energy 

Each time the bath and system are disconnected the potential energy of interaction is 
canceled; thus 

dEr _ _dEA ^dE^^dEA , . 

dt dt dt ~* dt 

where the last quantity represent the asymptotic value when S has reached its stationary 
state. Note that this depends on the initial interaction energy being zero (|6]) as well as 
on the conservation of energy during the coupling interval. The asymptotic interaction 
energy flow, ^§f-, is negative if the bath has more than one frequency and results from the 
coupling/decoupling process at each time step. As indicated, it corresponds to the work on 
the system S provided by an external source. 

4. Spin boson energy flows 

We mention corresponding results for the SB interaction. As for RW, ^f- — > exponen- 
tially. The flow of energy to the bath is 

dE A 



dt 



-7 2 5> fc | 2 u;(fc)<0. (B6) 



Once again, a constant flow of work must be provided by an external source during the 
coupling/decoupling process. 
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